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Abstract. The problem of construction of Barabanov norms for analysis of 
properties of the joint (generalized) spectral radius of matrix sets has been dis- 
cussed in a number of publications. In 18 , 21 the method of Barabanov norms 
was the key instrument in disproving the Lagarias-Wang Finiteness Conjecture. 
The related constructions were essentially based on the study of the geomet- 
rical properties of the unit balls of some specific Barabanov norms. In this 
context the situation when one fails to find among current publications any 
detailed analysis of the geometrical properties of the unit balls of Barabanov 
norms looks a bit paradoxical. Partially this is explained by the fact that Bara- 
banov norms are defined nonconstructively, by an implicit procedure. So, even 
in simplest cases it is very difficult to visualize the shape of their unit balls. 
The present work may be treated as the first step to make up this deficiency. 
In the paper an iteration procedure is considered that allows to build numer- 
ically Barabanov norms for the irreducible matrix sets and simultaneously to 
compute the joint spectral radius of these sets. 

1. Introduction. Let = {Ai, . . . , A^} be a set of real m x m matrices. As 
usual, for n > 1, let us denote by jz/" the set of all n-products of matrices from £/ ; 
s^'^ = /. For each n > 1, define the quantity 

Pn{£^) := max p{Ai^^ ■ ■ ■ A,;,AiJ, 

where maximum is taken over all possible products of n matrices from the set si ^ 
and p(-) denotes the spectral radius of a matrix, that is the maximal magnitude of 
its eigenvalues. Clearly, if n > r then some matrices in the product Ai^ ■ ■ ■ Ai^Ai^ 
will occur several times. The limit 

p(^) := limsup(p„(^))^^" 

n— >cxD 

is called the generalized spectral radius of the matrix set .s/ PlITTj . 

Similarly, given a norm j| • || in M'" , for each n> 1, define the quantity 

Pn{£/) ■■= max ••• Ai^AiJI, 
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where \\A\\, for a matrix A, is the matrix norm generated by the vector norm j| • | 
in R"\ that is = sup||^||^i \\Ax\\. Then the hmit 

pK) := limsup (p„(=^))^/" (1) 

does not depend on the choice of the norm || • || and is called the joint spectral radius 
of the matrix set £/ [33 . When r = 1 this definition coincides with the famous 
Gelfand formula for the spectral radius of a matrix [13 as in this case £/ ~ {A} is 
a singleton matrix set and /5„(j2/) — 

For matrix sets consisting of a finite amount of matrices, as is our case, the 
quantities p(^) and p{s!/) coincide with each other 5 and their common value is 
denoted as 

while the quantities /3„(^) and Pnis^) form lower and upper bounds, respectively, 
for the joint/generalized spectral radius: 

Pn{^) <P{^) = P{^) < Pn{-S^), V n > 0. 

This last formula may serve as a basis for a posteriori estimating the accuracy of 
computation of p{s^). The first algorithms of a kind in the context of control theory 
problems have been suggested in I6J, for linear inclusions in [2 , and for problems 
of wavelet theory in 8 -10 . Later the computational efficiency of these algorithms 
was essentially improved in 14,26 . Unfortunately, the common feature of all such 
algorithms is that they do not provide any bounds for the number of computational 
steps required to get desired accuracy of approximation of p{s!^) ■ 

Recently, in [19 explicit computable estimates for the rate of convergence of the 
quantities to p{A) end their extension to the case of the joint spectral 

radius were obtained. Probably, these results will help to make more constructive 
the problem of evaluating of p{s^) by the generalized Gelfand formula (1) . 

Some works suggest different formulas to compute p{s^). So, in [7^ it is shown 
that 



p{s^) ~ limsup max |tr(Ai^^ • • ■ Ai^Ai-^] 



where, as usual, tr(-) denotes the trace of a matrix. 
Given a norm || • || in M™, denote 

Kll :=niax||A||. 

Then the spectral radius of the matrix set =e/ can be defined by the equality 

pK)=infK||, (2) 



where infimum is taken over all norms in R [12(133 . and therefore 

for any norm || • || in R™. For irreducible matrix sets,/ the infimum in (2) is attained, 
and for such matrix sets there are norms || • || in R™, called extremal norms, for 
which 

K||<p(^). (3) 



matrix set .s/ is called irreducible if the matrices from is/ have no common invariant sub- 
spaces except {0} and . In 23 25 such a matrix set was called quasi-controllable. 
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In analysis of the joint spectral radius ideas suggested by N.E. Barabanov 
play an important role. These ideas have got further development in a variety of 
publications among which we would like to distinguish 35 . 

Theorem 1.1 (N.E. Barabanov). Let the matrix set = he irre- 

ducible. Then the quantity p is the joint (generalized) spectral radius of the set si 
iff there is a norm || • || in M™ such that 



Throughout the paper, a norm satisfying (4) will be called a Barabanov norm 
corresponding to the matrix set s/ . Note that Barabanov norms are not unique. 

Similarly, as is shown in [31 , Thm 3.3] and [32 , the value of p equals to p{£/) if 
and only if for some central-symmetric convex bodj-^ 5* the following equality holds 



where conv(-) denotes the convex hull of a set and pS := {px : x E S}. As is noted 
in [31J, the relation (5) was proved by A.N. Dranishnikov and S.V. Konyagin, so it 
is natural to call the central-symmetric set 5* the Dranishnikov-Konyagin-Protasov 
set. The set S can be treated as the unit ball of some norm || • || in M'' (recently 
this norm is usually called the Protasov norm). Note that Barabanov and Protasov 
norms are the extremal norms, that is they satisfy the inequality (3). In [28,29,36 
it is shown that Barabanov and Protasov norms are dual to each other. 

Remark that formulas (3), (4) and (5) define the joint or generalized spectral 
radius for a matrix set in an apparently computationally nonconstructive manner. 
In spite of that, namely such formulas underlie quite a number of theoretical con- 
structions (see, e.g., Ill|18l|21]|27l|35il36j) and algorithms [30 for computation of 



Different approaches for constructing Barabanov norms to analyze properties 
of the joint (generalized) spectral radius are discussed, e.g., in [15, 17 and 34, 
Sect. 6.6]. In T8,21^ the method of Barabanov norms was the key instrument 
in disproving the Lagarias-Wang Finiteness Conjecture. The related constructions 
were essentially based on the study of the geometrical properties of the unit balls 
of some specific Barabanov norms. In [16,17 the method of extremal polytope 
norms was the key tool in investigation of the finiteness properties of pairs of 2 x 2 
matrices. 

In this context the situation when one fails to find among current publications 
any detailed analysis of the geometrical properties of the unit balls of Barabanov 
norms looks a bit paradoxical. Partially this is explained by the fact that Barabanov 
norms are defined nonconstructively, by an implicit procedure. So, even in simplest 
cases it is very difficult to visualize the shape of their unit balls. The present work 
may be treated as the first step to make up this deficiency. 

In the paper, an iteration procedure is considered that allows to build numerically 
Barabanov norms for the irreducible matrix sets and simultaneously to compute the 
joint spectral radius of these sets. A similar iteration procedure is also discussed 
in [22 . 

The structure of the paper is as follows. In Introduction we give basic definitions 
and present the motivation of the work. In Section 2 the main iteration procedure is 



P\\x\ 



max II 



(4) 




(5) 



^The set is called body if it eontains at least one interior point. 
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introduced and Main Theorem stating convergence of this procedure is formulated. 
The iteration procedure under consideration is called the max-relaxation procedure 
since in it the next approximation to the Barabanov norm is constructed as the 
maximum of the current approximation and some auxiliary norm. In Section 3 proof 
Main Theorem is given. In Section 4, to build simplest examples, the max-relaxation 
scheme is adapted for computations with 2x2 matrices. Results of numerical tests 
are illustrated by two examples. At last, in concluding Section 5 we discuss some 
shortcomings of the proposed approach and formulate further unresolved problems. 

2. Max-relaxation iteration scheme. Given r,m > 1, let £/ ~ {Ai, . . . , Ar} be 
an irreducible set of m x m real matrices. 

Throughout the paper, a continuous function 7(t, s) defined for t,s > Q and 
satisfying 

7(t, t) = t, min{t, s} < j{t, s) < max{t, s} for t ^ s, 
will be called an averaging Junction. Examples for averaging functions are: 

/ \ i '\' s . I — , . 2is 
7(^, s) = l[t, s) = Vis, 7(i, s) = —— . 

Given some averaging function 7(-, ^ norm || • j|o in R™, and a vector e S M™ 
such that ||e||o = 1, construct recursively the norms || • ||„ and || • ||°, n ~ 1,2, . . ., 
in accordance with the following rules: 

MRl: if the norm \\ ■ ||„ has been already defined compute the quantities 

, maxi|jAix||„ _ . maxi||Aix||„ _ , 

MU ' '-^"^^ INU ' 7« = 7(P„,PJ; (6) 

MR2: define the norms \\ ■ and \\ ■ 

||a:||„+i = max |||.t||„, 7;^^ max || Aj2;||„| , (7) 
= ||.T||„+i/||e||„+i. (8) 

Remark that the number of operations needed to perform one step of algorithm 
MR1-MR2 is of order rm^v{£), where vi^e) is the number of operations needed to 
compute, for an arbitrary vector x G K™, the value of the norm \\x\\ with a relative 
accuracy e. In general, the value v{e) is of order e~™. So, the total number of 
operations needed to perform n steps of algorithm MR1-MR2 has the same rate 
of growth as nrm'^ e~"^ . 

Remark also, that the procedure (6) of calculation of resembles the technique 
of iterative approximation of the joint spectral radius suggested in 14 . 

Main Theorem. For any irreducible matrix set £/ , nonzero vector e G M™, initial 
norm \\ ■ ||o, and any averaging function 7(<, s), the sequences {p^} constructed by 
the iteration procedure MRl, MR2 converge to pi^si), and the sequence of norms 
II • II „ converges uniformly on each bounded set to some Barabanov norm || • ||* of the 
matrix set . Moreover, the sequence {p^} is nondecr easing, the sequence {pX\ is 
nonincreasing, and 

Pn < P{^) < Pt 

for all n ~ 1,2, ... , which provides an a posteriori estimate for the computational 
error of p{j^). 
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3. Proof of Main Theorem. Let us suppose that we managed to prove the fol- 
lowing assertions: 

Al: the sequences {p^} o,nd {p^} are convergent; 

A2: the limits of the sequences {p^} o,rid {p^} coincide: 

p = Urn p+ = Hm p-; 
A3: the norms \\ ■ ||° converge pointwise to a limit \\ ■ \\* . 

Then the function || • ||* will be a semi-norm in R"'. Moreover, by (8) each norm 
II • 11° meets the normalization condition ||e||° = 1. and hence 

||e|r = hm ||e||° =1, 

which implies ||a;||* ^ 0. Note also that due to (8) the norms || • ||° differ from || • ||„ 
only by numerical factors. Therefore, the quantities p^ can be defined as 

, maxi||Ajx||° _ . maxj ||Aia;||° 
p,T = max ;— r , p., = mm ;— r . (9) 

||x||° ^" ||.t||° ^ ' 

Then, passing to the limit in (9), one can conclude that the semi- norm ||a;||* satisfies 
the Barabanov condition 

pII^II* max ||ylia:||*. 

i 

But as shown in [21, Thm. 3], any semi-norm ||a:||* ^ satisfying the Barabanov 
condition for an irreducible matrix set is a Barabanov norm. 

Thus, under assumptions Al, A2 and A3, the iteration procedure (6)-(8) allows 
to build a Barabanov norm and to find the joint spectral radius of the matrix set 

So, to complete the proof of Main Theorem we need to justify assertions Al, 
A2 and A3 which will be done in Sections 3.1-3.6, In Section 3.1 we establish 
convergence of the sequence of norms {|| • ||°i} to some limit which allows to prove in 
Lemma 3.2 that Assertion A3 is a corollary of Assertions Al and A2. Section 3.2 
demonstrates that the quantities form the family of lower bounds for the 

joint spectral radius p of the matrix set while the quantities {p^} form the 
family of upper bounds for p. In Section 3.3 we prove that the sequences {p^} are 
bounded and monotone which implies the existence of the limits p~ = lim„_j.oo Pn 
and p"'" = lim„_i.oo Pn - \&s.t, in Sections 3.4 and 3.5 we prove that p^ = p+ which 
allows to justify in Section 3.6 the validity of Assertions Al and A2 and thus to 
finalize the proof of Main Theorem. 

3.1. Convergence of the sequence of norms {|| • ||°}. Given a pair of norms 
II • II' and II • II" in M™ define the quantities 

e-(IMIMMI")-min|^, e^(MIM|-||")=-xW;. (10) 

Since all norms in R™ are equivalent to each other, the quantities e^(|| • ||', || • ||") 
and e'^dl • ||', || • ||") are correctly defined and 

0<e-(||.||M|.||")<e+(||.|r,||.||")<oo. 

Therefore the quantity 

«'(iMiMi-in- ::;i|:||:;[i:||::; >i. dD 
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which is called the eccentricity of the norm || • ||' with respect to the norm || • |j" 
(see, e.g., 36 ), is also correctly defined. 

Lemma 3.1. Let j| • ||* &e a Barabanov norm for the matrix set £/ . Then 

ecc(||.||°,||.r) = ccc(||.||„,||.|r), Vn, (f2) 

and the sequence of the numbers ecc(|| • ilm il ' H*) nonincreasing. 

Proof. Note first that by (8) each norm || • ||° differs from the corresponding norm 
II • ||„ only by a nmiierical factor. From this, by the definition (fO), (ff ) of the 
eccentricity of one norm with respect to another, the equality (12) follows. 

Denote by p the joint spectral radius of the matrix set .s^/. Then, by definitions 
of the function e+(-) and of the Barabanov norm || • ||*, from (6), (7) we obtain: 

||a;||„+i = max|||a;||„, 7;;;^ max || Aj.t||„| < 

< e+(|| • ||„,|| • r)max{||a;|r, 7,7' max || A,.T|r } = 

= e+(|| • ||„, II • ||*)max{||x||*, -/^^^ p\\x\\*} . 

Therefore 

e+(M|„+i,||-|r)<e+(||.|U,||.|r)max{l, 7,7 V} ■ (13) 
Similarly, by definitions of the function e^(-) and of the Barabanov norm || • ||*, 
from (6), (7) we obtain: 

||a;||„+i = max|||a::||„, 7"^ max || A^a::||„| > 

> e"(|| • ||„,|| • ||*)max{||x|r, 7-1 max ||A,.t||*} = 

= e-(i|.|u,i|.r)max{iixir, 7„-viia;ir} • 

Therefore 

e-(||-|U+i,Mr)>e-(||.||„,||.|r)max{l, 7-V} • (W) 
By dividing termwisc the inequality (13) on (14) we get 

e+(i|.|u+i,i|.||*) ^ e+(i|.|u,i|.r) 

ecc( • „+i, • = _.|| II n— [TIT ^ -/|i II — riilT = ■ ■ )■ 

e (II ■ l|n+i,ll ■ Ir) e (II • ||„,|| • II*) 

Hence, the sequence {ecc(|j • ||„, || • ||*)} is nonincreasing. □ 
Denote by A'ioc(K™) the topological space of norms in M™ with the topology of 
uniform convergence on bounded subsets of R'" . 

Corollary 1. The sequence of norms {\\ ■ ||°} is compact in Niod^"^)- 

Proof. For each n and any x 7^ 0, by the definition (10) of the functions e'''(-) and 
e^(-) the following relations hold 



and then 



from which 



<|^<e+(||.||„,||.r). 



e-(M|„,||-r)<f^<e+(||.||„,||.|r), 



H|e||„<||x||„<ccc(||-||„,||-r)pi-||e||„. 



ecc • • r e r e 
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Since here by construction the norms {|| • ||°} satisfy the normahzation condition 
||e||° = 1, and by Lemma 3.1 ccc(|| • || • < ccc(|| • ||°, |1 • ||*), we have 

^ < 11-^11" <ccc(M|o, II- 



ecc(||-l|o,l|-|r) Ih 

Therefore the norms || • ||n, ?^ > 1, are equicontinuous and uniformly bounded on 
each bounded subset of K™. Moreover, their values are also uniformly separated 
from zero on each bounded subset of W" separated from zero. From here by the 
Arzela-Ascoli theorem the statement of the corollary follows. □ 

Corollary 2. // at least one of subsequences of norms from {\\ ■ ||°} converges in 
Afioc(K™) to some Barabanov norm then the whole sequence {\\ ■ ||°} also converges 
in 7Vioc(K'") to the same Barabanov norm. 

Proof. Let {|| • ||°^} be a subsequence of {|| • ||°} which converges in Afioc(]R™) to 
some Barabanov norm || • ||*. Then by definition of the eccentricity of one norm 
with respect to another 

ecc(|| • ||°^,|| • ID ^ 1 as/c^oo. 

Here by Lemma 3.1 the eccentricities ccc(|| • ||°, || • ||*) arc nonincrcasing in n, and 
then the following stronger relation holds 

ecc(||-||°,|l-|r)^l asn^oo. (15) 

Note now that by the definition (10), (11) of the eccentricity of one norm with 
respect to another 

— ^^^^<Si<ecc(||.||°,||.r), 



^^^Vll lln' II II J 11*^11 

from which by (15) it follows that the sequence of norms {j| • ||°} converges in space 
iVioc(R") to the norm || • ||*. □ 

Lemma 3.2. Assertion A3 is a corollary of Assertions Al and A2. 

Proof. By Corollary 1 the sequence of norms {|| • ||,°} has a subsequence {|| • ||,°^} 
that converges in space NiocQSJ") to some norm || • ||*. Then, passing to the limit 
in (9) as n = nfc — >■ oo, we get by Assertions Al and A2; 

max,; \\A^x\\* 

which means that || • ||* is a Barabanov norm for the matrix set £/. This and 
Corollary 2 then imply that the sequence {|| • ||°} converges in space A'^ioc(K™) to 
the Barabanov norm || • ||*. Assertion A3 is proved. □ 
So, in view of Lemma 3.2 to prove that the iteration procedure (6)-(8) is con- 
vergent it suffices to verify only that Assertions Al and A2 hold. 

3.2. Relations between and p. The following lemma provides a way to esti- 
mate the spectral radius of a matrix set. 

Lemma 3.3. Let a, (3 be numbers such that in some norm \\ ■ \\ the inequalities 

a\\x\\ < max ||^i.'j;|| < /3||a;||, 

hold. Then a < p < P, where p is the joint spectral radius of the matrix set jz/. 
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Proof. Let || • ||* be some Barabanov norm for the matrix set Since all norms 
in R™ are equivalent, there arc constants a~ > and tr"'" < oo such that 

(T^llxll* < ||a;|| < cr+|lx||*. (16) 

Consider, for each fc = 1,2,..., the functions 

Afc(a;) = max \\A,^ . . . A^^Ai^x\\. 

I<ii,i2;---,ik<r 

Then, as is easy to see, 

a'^llxll < Afc(.T) </3^||a;||. (17) 
Similarly consider, for each fc = 1, 2, . . ., the functions 

Afe(a;)= max \\A,...A,A,x\\*. 

For these functions, by definition of Barabanov norms the following identity hold 

A^(x)^/||x|r, (18) 

which is stronger than (17). 

Now, note that (16) and the definition of the functions Afc(x) and A^(a;) imply 

a-AUx) < Ak{x) < a+AUx). 

Then, by (17), (18), 

cr+ cr 

from which the required estimates a < p < (3 follow. □ 
So, Lemma 3.3 and the definition (6) of imply that the quantities {p~} form 
the family of lower bounds for the joint spectral radius p of the matrix set £/, while 
the quantities {Pn} form the family of upper bounds for p. This allows to estimate 
a posteriori errors of computation of the joint spectral radius with the help of the 
iteration procedure (6)-(8). 

3.3. Convergence of the sequences {p^}- Estimate the value of max^ 
By definition, 

max ll^jxlln+i = max < max < || Aia;||„, max || Ai^jx||„ 

= max < max ||^ia:||„, 7,7^ maxmax ||^iy4ja;||„ 

[ i 3 i 

Here by the definition (6) of the quantities p^ the right-hand part of the chain of 
equalities can be estimated as follows: 

p„ max|||a;||„, 7^^ max Pjx||„| < 

< max < max Ij^ixjln, 7""'^ maxmax ||^iylja:;j|„ > < 
I i ' 3 I J 

p+ max |||a::||„, 7,7^ max || AjxH 

Therefore, by definition of the norm j|.xj|„+i, 

P,;lla:^||n+i < ma.x\\Aix\\n+i < p+|lx||„+i, 



ITERATIVE BUILDING OF BARABANOV NORMS 



9 



from which 

_ max, P,a;||„+i + w /n 

and then, 

Pn < Pn+l ^ Pn+l < Pn ■ 

So, the fohowing lemma holds. 

Lemma 3.4. The sequence {p^} is bounded from above by each member of the 
sequence {p^} and is nondecreasing. The sequence {p^} is bounded from below by 
each member of the sequence {p~} and is nonincreasing. 



In view of Lemma 3.4 there are the limits 

p^ = lim p+ = lim p+, 7 lim 7„ = lim 7(p^,p+), 



n— >oo 

where 

< 7 < P^, 

which means that Assertion Al holds. Hence, to prove that the iteration procedure 
(6)-(8) is convergent it remains only to justify Assertion A2: = p+. 

To prove that p" = p^, below it will be supposed the contrary, which will lead 
us to a contradiction. 

3.4. Transition to a new sequence of norms. To simplify further constructions 
we will switch over to a new sequence of norms for which the quantities will be 
independent of n. 

By Corollary 1 the sequence of the norms || • ||° is compact in space A'loclK™). 
Hence, there is a subsequence of indices {rifc} such that the the norms || • ||°^ = 
II ■ ll"fc/l|s|l"fc converge to some norm || • ||5 satisfying the normalization condition 
Hell* = 1. Then, passing to the limit in (9), by Lemma 3.4 we obtain: 

, maxi ||A,a:||;; _ . max, ||A,;a-||J _ , 

P - II.II5 ' P = \\A\l ' ' ' 

Now by induction the following statement can be easily proved. 

Lemma 3.5. For each n ~ 0,1,2,..., the sequence of the norms \\ ■ ||nfc+n/||e||„j. 
converges to some norm \\ ■ ||*. Moreover, for each n ~ 0,1,2,..., we have the 
equalities 

max, ||A,a:||* ^ . max, ||A,a:||' _ 

max — = p , mm t—t. —p , (19) 

x/O ||x||' .#0 ||x||' ' ^ ' 

and the recurrent relations 

||x||:+i=max{||a;||:, 7"! max || A,x||:} . (20) 

3.5. Sets a;„. Define, for each n = 0, 1, 2, . . ., the set 

w„ = {:e e M'" : p-||a;||,* = max || A,x||,* } . (21) 

By (19) bJn is the set on which the quantity 

maxj ll^jxll* 

attains its minimum. 

Lemma 3.6. If x £ w„ then ||a;||*_|_]^ = ll^^ll' - 
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Proof. The statement of the lemma is obvious for x = 0. So, suppose that x € uJn, 
X ^ Q. In this case (21) and the inequahtics < p+ imply 

max||A,a;||: = p-i|a;||:<7lN||: 

or, what is the same. 

From here by the definition (20) of the norm || • W^+i wc get the required equality: 

=max|||x||', max |j Ax||;| = ||a;||*. 

The lemma is proved. □ 

Lemma 3.7. // < p+ then uJn+i ^ for each n = 0, 1, 2, . . .. 

Proof. Let x € Wn+i- If a; = then clearly x S w„, so suppose that x ^ Q. By 
definitions of the set a;„+i and of the norm || • ||* the following equalities take place: 



max II Ai.T||*^j = max |max ||| Aixjl' , 7 ^ max ||AjAi.T||* 

= max < max || A,ia:||' , max || AjAix||* > = 
t ' J 

= P"l|a^ll'+i = P"max|||a:||', 7"^ max || A,x||' | . (22) 

Let here 

||x||: <7"'max||A,a;||,'. (23) 
Then from (22) it follows that 

max<^ max ||A,.t||', 7"^ max || A^ Aia;||' \ =p~||2^||' +1 = 7" V~ max || A,(.t||* . 

1^ ^ «j' J ^ 

But by the conditions of the lemma p~ < . Then 7 ~ 7(p^,p+) > p^, and the 
right-hand part of the above equalities is strictly less than max^ ||AiX'||' . A contra- 
diction, since the left-hand part of the same equalities is no less than max, ||Aix||* . 

The above contradiction is caused by the assumption (23), and therefore it is 
proved that the condition a; 7^ G ujn+i implies the strict inequality 

> 7"^max||A,a:||*. 
In this case from ( 22 ) it follows that 



max<^max||A,x||', 7 ^ max || Aj A,x||' ^ = p ||a;||*. (24) 

Let us show that the equality (24) implies 

max||A,x||: -p-||a;||,'. (25) 

Indeed, supposing the contrary, by definition of the quantity p~, there should be 
vahd the strict inequality max^ ||Aia;||* > p~||x||'. Then the left-hand part of the 
equahty (24) should be strictly greater than p~||a;||*, that is greater than the right- 
hand part of the same equality, which is impossible. This last contradiction shows 
that the equality (25) holds as soon as a; ^ G Wn-i-i) which means by (19) that 

X £ bJn- □ 
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Corollary 3. If p < /?+ then uj := nri>ow„ 7^ and 

ll-^llo = l|.^llI = --- = l|2:||' Vx^Oeuj. (26) 

Proof. By Lemma 3.7 {a;„} is the family of embedded closed conic sets. Then the 
intersection uj of these sets is also a closed conic set such that lj ^ {0}. 

By definition of the set w, if x E uj then x e for every integer n > 0. Hence, 
by Lemma 3.6 |la::||*+i — ll^;]!*, from which the equalities (26) follow. □ 

3.6. Completion of the proof of Assertion A2. By Corollary 3 there is a non- 
zero vector g on which all the norms || • ||,* take the same values: 

ll.9llo = ll.9llI = --— Il5li: = ---- 

Then, due to uniform boundedness of the eccentricities of the norms |j • ||* with 
respect to some Barabanov norm || • || * (this fact can be proved by verbatim repetition 
of the analogous proof for the norms || • ||„), the norms || • ||* form a family which is 
uniformly bounded and equicontinuous with respect to the Barabanov norm || • ||*: 

3M±e(0,cx)): ti-\\x\\* <\\x\\:^<fi+\\x\\*, n = 0,l,2,.... 

Hence, by the Arzcla-Ascoli theorem the family of norms {|| • ||*} is compact in 
iVioc(K™). 

From the definition (20) of the norms || • ||,* it follows also that 



7 1 max||A,x||*| 



> lla^ll!. 



Then the norms || • ||* are monotone increasing in n and bounded (with respect to 
the Barabanov norm || • ||*) and therefore they pointwise converge to some norm 
II • II*. Moreover, since the family of norms {|| • ||*} is equicontinuous with respect 
to the Barabanov norm || • ||*, the norms || • ||* converge to the norm || • ||* in space 

^loc(M™). 

Now, passing to the limit in the relations 

Ii2;||*+i =max|||x||', 7"^max||Aix||'| > 7"^ max || Aia::||* , 
which follow from (20), we obtain 

||a;|r > 7^^max|lAja;||*. 

i 

From here 

max," IjAixll* 

max , "„ " < 7. 27 

^/o ||x||' - ' ' 

On the other hand, passing to the limit in the first relation of (19), we obtain 

max,- ||A,a;||* , 
max , "„ " =p+. 28 

Relations (27) and (28) imply the inequality < j which contradicts the as- 
sumption < p"*" because by definition of the function 7(-, •) the condition p~ < p"*" 
implies the inequality 7 — 7(p~,/9+) < p+. 

The obtained contradiction completes the proof of the equality p^ = as well 
as of the convergence of the iteration procedure (6)-(8). 



set X is called conic if together with each its point x it contains the ray {tx : t > 0}. 
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4. Computational scheme for 2x2 matrices. Let = {^i, . . . , Ar} be a set 

of real 2x2 matrices 



(i) Ji) 



11 


ai2 


(i) 
21 


(i) 
0-22 



A, = 



Let (r, (f) be the polar coordinates in . Then, for a vector x £ with 
Cartesian coordinates x = {a;i,a;2}, we have 

X = {r cos iy9, r sin ip} 

and 

r = r(a;) — \J x\ + x\^ Lp = ^{x) ~ arctan (2;2/a;i) . 

Define, for an arbitrary norm || • ||, the function 

R{lp) = ||{cos.^,sin(p}||. 

Then the norm ||a;|| of the vector x with polar coordinates {r,Lp) is determined by 
the equality 

\\x\\=rR{ip), (29) 

and the unit sphere in the norm || • || is determined as the geometrical locus of the 
vectors x polar coordinates of which satisfy the relations 

1 



rR{tp) — 1 or r = 



R{^y 



see Fiff. 1 





f x = {r, tp) 




^J^\x\\ = rR{ip) 




/ l/R{pj\ 


= 1 <^ rR{ip) 


= 1 ] 



Figure 1. Definition of the function R{ip). 



Now, let Rn{f) be the function defining in the polar coordinates the graph of 
the unit sphere |jx|j„ = 1 of the norm || • ||„ determined by the iteration procedure 
(6j-(8). Rewrite the relations (6)-(8) in terms of the functions Rn{<p)- To do it 
we should express the quantities ||j4ia;||„, i — 0, 2, . . . , r, in terms of the functions 

By (29) 

\\Aix\\n = r{A,x)Rn{^iAix)). 
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Here by definition of the matrix Ai 

r{Aix) = rHi{ip), 

where 



rii[ip) ~ y [<^ii cosifi + a\2 s,inip I + \a2i cos + smip I 



Similarly, by definition of the matrix Ai 
where 



^i{^p) = arctan 



cos tp + sm ip 
^ aj^/ cos (p + a]^2 sm ip ^ 



From the obtained relations it follows that the first two equalities in (6) take the 
form 

+ i7,((p)i?„($,((^)) _ . i7,(^)i?„($,((p)) 
pj^ ^ max max — — , p„ = mm max — — , 

or, what is the same, 

p+=max— — , p„ = mm (30) 
V RnW) V n,i[(p) 

where 

KCV') = max{i/,(^)i?„($,((/p))} . (31) 

i 

The relations (7) take the form 

rRn+i{tp) = max{rR„{'p), n-^R*„{'p)} 

or, equivalently, 

Rn+iiv) = max {i?„(^), 7,7'K(V')} (32) 
and the normalization condition (8) takes the form 

reRn+l((Pe) 

where (re, ipe) are polar coordinates of the vector e. Taking in place of e the vector 
with polar coordinates (1,0) the normalization condition can be rewritten in the 
form 

i?:+iM = §^. (33) 

So, the max-relaxation iteration scheme can be represented as follows. Given 
an averaging function 7(-, •), set i?o(<p) = li and build recursively the 27r-periodic 
functions i?„(v?) and i?° ((/?), n = 1, 2, . . ., in accordance with the following rules: 

MRl: regarding the function Rn{'p) already known compute the numerical val- 
ues and p^ by formulas (30), (31) and set 7„ = 7(p^, pt); 

MR2: define Rn+iiip) by (32) and R°^_^i{(p) by (33), and then determine the 
norm \\ ■ W^+i as II^^IIJ^+i = ri?°_|„]^ (</?), where (r,ip) are the polar coordinates 
of the vector x. 
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2 
1.5 
1 

0.5 


-0.5 
-1 
-1.5 

-2 






= P 


- - Uix\\ 


= P 


Ps^ll 


= P 




= 1 




Figure 2. Examples of computation of Barabanov norms for 2x2 matrices. 



Example 1. Consider the family 



{^1, A2} of 2 X 2 matrices 



The functions 3>i (lys) , Hi (1^9) , i?„ {^p) , i?* (1^9) were chosen to be piecewise linear with 
3000 nodes uniformly distributed over the interval [— 7r,7r]. It was needed 13 iter- 
ations of algorithm MR1-MR2 with the averaging function 7(t, s) = imple- 
mented in MATLAB to compute the joint spectral radius p(^) with the absolute 
accuracy lO"''. The computed value of the joint spectral radius is p(^) — 1.617. 
The computed unit sphere of the Barabanov norm || • ||* after the 13th iteration of 
algorithm MR1-MR2 is shown on Fig. 2 on the left. 

As is seen from Fig. 2, in Example 1 the sets = p and ||A2a:;|| = p have 

exactly 4 intersection points. This was theoretically proved in 18,21 for the case 
when one of the matrices Ai , A2 is lower triangle and the other is upper triangle, 
and their entries are nonnegative. In 18,21 this fact was one of key points in 
disproving the Finiteness Conjecture. We do not know whether this fact is valid 
in a general case or not, but numerical tests based on algorithm MR1-MR2 with 
several dozens pairs of matrices Ai,A2 testify for this fact. 



Example 2. Consider the family 



= {Ai, A2, A3} of 2 X 2 matrices 

A^ = 



0.8 
-0.6 



0.6 
0.8 



1 

-0.4 





1.3 



Here the functions ^'i(</3), Hi{(p), Rn{(p), Rnif) were also chosen to be piecewise lin- 
ear with 3000 nodes uniformly distributed over the interval [— 7r,7r]. It was needed 
31 iterations of algorithm MR1-MR2 with the averaging function jit, s) = im- 
plemented in MATLAB to compute the joint spectral radius p{£/) with the absolute 
accuracy 10~^. The computed value of the joint spectral radius is p{£^) = 1.347. 
The computed unit sphere of the Barabanov norm || • ||* after the 3 Id iteration of 
algorithm MR1-MR2 is shown on Fig. 2 on the right. The MATLAB code for this 
Example can be found in [20 . 



ITERATIVE BUILDING OF BARABANOV NORMS 



15 



5. Concluding remarks. The max-relaxation algorithm suggested in the paper 
allows to calculate the joint spectral radius of a finite matrix family (of arbitrary 
matrix size and arbitrary amount of matrices in the set) with any required accuracy 
and to evaluate a posteriori the computational error. Still, this algorithm gives rise 
to a set of open problems. 

Problem. While the quantities {p^} provide a posteriori bounds for the accuracy 
of approximation of p{£^) the question about the accuracy of approximation of the 
Barabanov norm \\ ■ \\* by the norms \\ ■ ||° is open. 

It seems, the difficulty in resolving this problem is caused by the fact that in gen- 
eral the Barabanov norms for a matrix family are determined ambiguously. Namely 
to overcome this difficulty we preferred to consider relaxation algorithm instead of 
direct one. Moreover, as can be shown, both theoretically and numerically, if to set 
||a;||„+i = maxi in (7) then the obtained direct computational analog of 

algorithm MR1-MR2 may turn out to be non-convergent. 

Problem. The question about the rate of convergence of the sequences {Pn} o,f^d 
to the joint spectral radius is also open. 

Remark also that in this paper mainly the algorithm for building of Barabanov 
norms rather than its computational details was studied. The numerical aspects of 
implementation of this algorithm require additional analysis. 

Problem. An estimation of the computational cost of the max-relaxation algorithm 
is required. The dependance of the algorithm on parameters r, m and the choice of 
the averaging function j(t, s) is acute, too. 
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